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Abstract- The present study deals with the constant and 
variable stiffness of the human leg joints, during the normal 
walking. A seven-link segment model was employed to represent 
the anterior-posterior motion study of the human body 
dynamics. The segments were interconnected by damped 
torsional springs for modeling the joint muscle-tendon complex. 
Lagrangian mechanics approach employed to representing 
equations of motion. The stiffness and damping for ankle, knee 
and hip joints can be estimated using numerical solution and 
optimization procedure. Damping was found negligible. There 
was closer agreement between constant joint stiffness given in 
the literature and the investigated values for stiffness 
coefficients from proposed model. The model presented may be 
provide an effective tool for future designing an active 
exoskeleton system to assist a paraplegic subject to walk by 
controlling the additional hip, knee and ankle joint stiffness and 
damping. 

I. Introduction 

Researches into human gait have a wide range of application 
in medicine, ergonomics, sport science and technology. Most 
often methods of multibody dynamics are used when 
investigation is focused on mechanical aspects of the gait [1- 
3]. Some multibody models for gait simulation represent 
bipeds of a very simple, not anthropometric kinematic 
structure [4], other models are limited to two dimensions [5- 
8]. Walking is a special kind of the multibody system motion 
with the occurrence of impacts, friction and slipping. These 
entire phenomenon observed during the foot-ground contact 
have to be introduced into the model to make it realistic. In 
some multibody models of the biped the foot-ground contact 
is modeled by the additional kinematic joints [5]. 

With the development of the biomechanical models of 
human body motion, it has become possible to simulate 
human performance in order to gain insight into 
intermuscular coordination and to elucidate control strategies 
of the musculoskeletal systems. A common method to deal 
with this type of problems is to lump together elements of the 
human body, e.g., muscles, tendons, ligaments, bones and 
joints so that the overall musculoskeletal system is 
represented as a damped elastic mechanism. 

Several methods describing the landing phase of running, 
hopping, or jumping can be found in the literature [9-13]. 
These models are usually characterized by the presence of 
elastic springs and viscous dampers, with constant properties 
and provide a reasonable prediction of the maximal vertical 
foot-ground reaction force. Depending on the range of joint 
flexion and on the frequency of motion, a considerable 
amount of elastic energy can be stored and reused. It has been 
shown that the dissipated energy in muscles increase when 
the amplitudes of joint movements are increased [14]. During 
the ground contact period of running, hopping and trotting, 
the leg was modeled as a single linear spring or in terms of 



the leg joints as constant torsional springs for the ankle, knee 
and hip joints, with no damping. The stiffness of these 
springs are considered the average value as the ratio between 
overall force or moment change to overall vertical 
displacement or angle change for the leg and joint stiffness, 
respectively [12, 15-17]. Physiologically, the conception of 
constant mechanical stiffness may not be applicable. For 
instance, muscular activation, behaved to be directly related 
to joint stiffness, varies during the stance phase in gait. 
Human joints and muscle-tendon complex have the ability to 
store, absorb and release the energy; therefore its stiffness 
generally depends on the activation level of the muscle. Thus 
it can be expected that the joint stiffness and damping is 
nonlinear in nature and that a model accounting for these 
facts may improve the human dynamical modeling and 
performance evaluation. 

Recently, there has been an increasing interest in the 
ability to control paraplegic stance. If individuals with 
paraplegia were able to stand comfortably for a prolonged 
period of time, it would have many therapeutic, psychological 
and practical advantages. Much of the researches regarding 
the control of paraplegic standing have been focused on the 
use of artificial stimulation to active paralyzed muscles and to 
partially restore motor function [18]. Activation of muscles 
either around the ankle, or the hip, or both assists the 
paralyzed individual in keeping the body's center of gravity 
within its base of support when stability is perturbed [19]. 
Paraplegic standing supported by FES-controlled ankle 
stiffriess is investigated in [20-22], in this study the knee and 
hip of the paraplegic subject were held in extended positions 
by a mechanical apparatus, and the closed loop control 
system maintain the appropriate level of stiffness around the 
ankle joint. In the other study, evaluation of the additional hip 
joint stiffriess on crutch supported paraplegic stance is 
investigated [23-24]. In a study illustrated in [25], the role of 
combined hip and ankle strategy during standing was 
implemented in a model as an optimal feedback postural 
control system. 

The main goal of the present work was to study the 
stiffriess and damping characteristics of the human joints in 
normal walking. The present study was thus pretended to 
provide an insight to realistic joint characteristic using the 
experimental gait data. It may be provide an effective tool for 
future designing an active exoskeleton system to assist a 
paraplegic subject to walk by controlling the additional hip, 
knee and ankle joint stiffness and damping. 

II. Mathematical model of the human body in walking 

A seven-link segment model is used to represent the 
human body gait dynamics as depicted in Fig. 1(a). The seven 
segments represent the left and right thighs, the shanks and 
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the feet, in addition, the head, arms and trunk (HAT) are 
illustrated by one segment. The mass, position of center of 
mass, moment of inertia, and length of each of the segments 
are scaled to the dimensions of a human (body mass 80 kg 
and 1.80 m height) using anthropometric data given in [26]. 
The segments are connected by torsional spring and damper 
for modeling the joint muscle-tendon complex. The model is 
planar and accounts only for anterior-posterior motion study 
of the body. For convenience, it is assumed that, initially, the 
right foot is forward with the right heel striking the floor. The 
joint moments are denoted by subscripted M's, with the first 
subscript indicating the right or left side of the body, and the 
second subscript indicating the ankle, knee and hip joints. 
Each of the feet is represented by a triangle of appropriate 
shape, as depicted in Fig. 1(b). 

The different phases of walking require two basically 
different mathematical representations. The different phases 
of walking are [3, 27]: 

Phase 1: From right-heel- strike (RHS) to left-toe-off 
(LTO). Here the right heel and left toe are pivoting on the 
floor. 

Phase 2: From LTO to right-foot-flat (RTF). Only the 
right heel is pivoting on the floor. 

Phase 3: During RFF. Here, the right foot is not moving, 
and the remaining six segments are pivoting on the right 
ankle. 

Phase 4: From RHO to left-heel- strike (LHS). Only the 
right toe is pivoting on the floor. 

During Phase 1, where six variables describing the two 
lower limbs are not independent, therefore as suggested in 
[27], the spring-damper combination in the leg joints plays 
vital role in the dynamic modeling, but during the swing 
phase. Phase 2, 3 and 4, the role of these elements were 
insignificant. In the present study the equations of motion for 
Phase 1 were considered and the detailed equations of motion 
for the latter cases may be found in [27]. 



dt 



'dL^ 
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dL ^ 
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/ = 1,2,...,7 



(2) 




The detailed formulation of the Lagrangian in Equation (1) 
and (2) are omitted because the intermediated expressions are 
very long, instead only the final equation of motion will be 
given. 

The equations of motion for double support case, RHS to 
LTO, arranged such that only the second derivatives are kept 
in the left hand side, has the following form: 



[«,]^=Z^ / = 1,2,...,7;7 = 1,2,...,7 

(2y =c^cos(6>,-6>^) 

7 

h=-Yc,a':^m(a^-a^+R 

7=1 



(3) 
(4) 



Since the Cjj are symmetrical, this makes the a,y symmetrical 
also. On the right hand side of the Equation (3), the bj are 
contain moments caused by the spring-damper combination, 
moments caused by gravity, the applied joint moments and 
the moments caused by centrifugal forces acting on the limbs. 
The expressions for Cjj and bj are given in Table L Also, the Cjj 
not given in the Table I are zero and all are symmetrical. 

Tablet 

c„ and R, COEFFICIENTS IN EQUATIONS (3) AND (4) 



Fig. 1 : Link segment model (a) limb angles and limb 
length variables (b) details of the foot variables 

Lagrangian mechanics provides the most direct approach" 
for obtaining equations of motion in this case. The 
Lagrangian mechanics approach requires the formulation of 
the Lagrangian given by: 

L = T-V (1) 

Lagrangian written in terms of the independent variables, 6^, 
therefore the resulting equation are given as follows: 
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in. Estimation of the joint stiffness and damping using 
Numerical solution 

In normal gait there is a natural symmetry between the left 
and right sides, this is graphically illustrated in Fig. 2. 
Therefore RHS to LTO takes place during 50% to 62% of 
gait cycle. Dynamic simulation of the proposed model 
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requires the knowledge of human kinetics and kinematics 
data. Joint angles and torques for a typical walking cycle 
were obtained in the form of independently collected Clinical 
Gait Analysis (CGA) data. CGA angle data is typically 
collected via human video motion capture. CGA torque data 
is calculated by estimating limb masses and inertias and 
applying dynamic equations to the motion data. Given the 
variations in individual gait and measuring methods, two 
independent sources of CGA data [31-32] were utilized for 
the analysis and estimation of the leg joint stiffness and 
damping. The anthropometric data required for simulation of 
the Equation (3) are specified in Table II [27]. In this case 
spring-damper coefficients of the leg joints are unknown in 
the proposed model. Numerical solution of the equations of 
motion in RHS to LTO phase result stiffness-gait cycle 
percent curves for ankle, knee and hip joints. Simulation 
results demonstrate angle dependent stiffness properties for 
joints during normal walking. Ankle, knee and hip stiffness 
profiles are depicted in Figs. 3 to 5 respectively. Damping 
coefficients of the joints were found negligible. Average 
values of investigated leg joints stiffiiess using numerical 
solution method are given in Table III and were compared 
with average values given in the literature [28-30]. 
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Fig. 2: The time spend on each limb during the normal gait cycle [3] 

Table ii 

ANTROPOMETRIC DATA FOR HUMAN BODY SEGMENT 



A 


Ankle to CG distance 


0.116 


B 


Heel to ankle distance 


0.140 


C 


Heel to CG distance 


0.189 


D 


Ankle to shank CG distance 


0.247 


E 


Shank CG to knee distance 


0.188 


F 


Knee to thigh CG distance 


0.227 


G 


Thigh CG to hip distance 


0.173 


H 


Hip to CG of head, arms and trunk (HAT) distance 


0.325 


a 


Angle between sole and ankle measured at the toe 


2.621 


p 


Angle between sole and ankle measured at the heel 


0.967 


r 


Angle between sole and CG measured at the heel 


0.310 


m, 


Foot mass 


1.159 


If 


Foot moment of inertia about CG 


0.010 


Ms 


Shank mass 


3.717 


Is 


Shank moment of inertia about CG 


0.064 


mt 


Thigh mass 


7.994 


It 


Thigh moment of inertia about CG 


0.133 


rrih 


Hat mass 


54.2 


h 


HAT moment of inertia about CG 


3.591 


w,„ 


Total body mass 


79.9 


^ 


Gravitational constant 


9.807 







_ ^ _ with CGA data given in [31] 
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Fig. 3: Ankle stiffness-Gait cycle percent during normal walking 
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Fig. 5: Hip stiffness-Gait cycle percent during normal walking 

IV. Estimation of the joint stiffness and damping using 

OPTIMIZATION PROCEDURE 

The main goal in this section is estimation of the joints 
stiffness and damping using optimization procedure. Similar 
to previous section two independent sources of CGA data 
[3 1-32] were utilized for the analysis. 
Elastic and damping torques are obtained by: 



[yM„=Pc,dcu, 



J = l,...,n , / = 1,2,3 



(5) 



here index / indicated the joint and n is the length of the joint 
moment vector. The stiffness and damping coefficients in 
Equation (3) can be solved by parameter estimation using 
optimization procedure. Minimization of the following 
objective function Jj was thus performed for each of the joints 
(ankle, knee and hip). 



7=1 



-M,y 



i = 1,2,3 



(6) 



here T- is CGA torque data. Four inequality constraints were 

applied as follows. The first two assign positive values for 
stiffness and damping, thus: 



C>0 



a, 2,3 



(7) 



The two next constraints were related to the potential elastic 
energy Est of the spring, and dissipating energy E^j of the 
damper, thus we have: 






j = \,...,n , z = 1,2,3 



(8) 



Parameter identification was performed by using Quadratic 
Programming (QP) method [33] to estimate the stiffness and 
damping of the ankle, knee and hip joints. Damping 
coefficients of the joints were found negligible using this 
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method too. Investigated values for joint stiffness coefficients 
using optimization procedure are given in the Table III. As 
illustrated in this table, there was closer agreement between 
constant joint stiffness given in the literature and the 
investigated values for stiffness and damping coefficients 
from proposed model using two different methods. 

Table hi 

COMPARISON OF THE LEG JOINT STIFFNESS VALUES IN 
THIS STUDY AND VALUES GFVEN IN THE LITRATURES 





Ankle 

(Nm/°) 


Knee 

(Nm/°) 


Hip 

(Nm/°) 


Stiffness given in [28] 


4.69 


2.67 


2.23 


Stiffness given in [29] 


4.23 


5.65 


3.96 


Stiffness given in [30] 


6.13 


5.03 


6.23 


Optimization results with CGA data [31] 


3.11 


2.12 


3.60 


Optimization results with CGA data [32] 


4.87 


3.39 


2.54 


Numerical solution with CGA data [31] 


3.32 


3.01 


5.62 


Numerical solution with CGA data [32] 


4.2 


2.50 


5.34 



IV. Conclusions and discussions 

The basic outcome of this study is estimation of leg 
joints stiffness and damping during normal walking. The 
value of the joints' damping coefficients was very small and 
negligible. There was closer agreement between constant 
joint stiffness given in the literature and the investigated 
values for stiffness and damping coefficients from proposed 
model using two different methods. FES-compensated ankle, 
knee and hip stiffness has been utilized to stand up paraplegic 
patients. By investigating the joints' spring-damper 
coefficients during normal walking, the presented model 
provides an effective tool for future designing of active 
artificial legs, orthosis or development of new walking 
control strategies with exoskeleton systems for paraplegic 
patients. 

References 

[I] Winter D.A., "Biomechanics of Human movement", John Wiley & Sons 
Inc., 1979. 

[2] Winter D.A., "Biomechanics and Motor Control of Human Movemenf , 
Second Edition, John Wiley & Sons Inc., 1990. 

[3] Vaughan C.L., Davis B.L., O'Conner J., "Dynamics of Human Gaif , 
Human Kinetics Publisher, Champaign, Illinois, 1992. 
[4] Bourassa P., Prince F., Meier M., Roy Y., "A Simple Walking Model; 
Simulation and Animation", Proc. of 6^^ International Symposium on 
Computer Simulation in Biomechanics, Tokyo, 1997, pp. 1-4. 
[5] Pandy M.G., Berme N., "Synthesizes of Human Walking: a Planar model 
for Single Support", Journal of Biomechanics, Vol. 21, 1988, pp. 1053-1060. 
[6] Koopman B., Grootenboer H.J., de Jongh H., "An Inverse Dynamics 
Model for the Analysis, Reconstruction and Prediction of Bipedal Walking", 
Journal of Biomechanics, Vol. 28, pp. 1369-1376, 1995. 
[7] Ouezdou F.B., Bruneau O., Guinot J.C, "Dynamic Analysis Tool for 
Legged Robots", Multibody System Dynamics, Vol. 2, No.l, 1998, pp. 369- 
391. 

[8] Wojtyra M., "Application of Human Gait Patterns to Biped Gait Control- 
Computer Simulation", Proc. of the 12^^ CISM-IFToMM Symposium on 
Theory and Practice of Robots and Manipulators, pp. 393-400, CISM, Udine, 
1998. 

[9] Mizrahi J., Suzak Z., "Elastic and Damping Response of the Human Leg 
to in Vivo Impact Forces", ASME Journal of Biomech. Eng. Vol. 104, 1982, 
pp. 63-66. 

[10] Ozguven H., Berme N., "An Experimental and Analytical Study of 
Impact Forces During Human Jumping", Journal of Biomechanics, Vol. 21, 
1988, pp. 1061-1066. 

[II] Kim W., Voloshin A.S., Johnson S.H., "Modeling of Heel Strike 
Transient during Running", Journal of Human Movement Science, vol. 13, 
1994, pp. 221-244. 



[12] Farley C.T., Houdijik H.H.P., Van Strien C, Lourie M., "Mechanism of 

Leg Stiffness Adjustment for Hopping on Surface of Different Stiffness", 

Journal of Applied Physiology, Vol. 85, 1998, pp. 1044-1055. 

[13] Meghdari A., Aryanpour M., "Dynamic Modeling and Analysis of the 

Human Jumping Process", Journal of Intelligent and Robotics Systems, Vol. 

37, 2003, pp. 97-115. 

[14] Bosco C, Komi P.V., "Mechanical Characteristics and Fiber 

Composition of Human Leg Extensor Muscles", Eur. Journal of Applied 

Physiology, Vol. 41, 1979, pp. 275-284. 

[15] Farley C.T., Morgenorth D.C., "Leg Stiffness Primarily Depends on 

Ankle Stiffness During Human Hopping", Journal of Biomechanics, Vol. 32, 

1999, pp. 267-273. 

[16] Johnson C, Hull M.L., "Parameter Identification on the Human Lower 

Limb under Dynamic, Transient Torsional Loading", Journal of 

Biomechanics, Vol. 21, No. 5, 1988, pp. 401-415. 

[17] Mote CD., Lee C.W., "Identification of Human Lower Extremity 

Dynamics in Torsion", Journal of Biomechanics, Vol. 15, No. 3, 1982, pp. 

211-222. 

[18] Kralj A., Bajd T., "Functional Electrical Stimulation: Standing and 

Walking after Spinal Cord Injury", Boca Raton, FL: CRC 1989. 

[19] Nashner L.M., McCollum G., "The Organization of Human Postural 

Movements: a Formal Basis and Experimental Synthesis", Behav. Brain Sci., 

Vol. 8, 1985, pp. 135-172. 

[20] Jaime R.P., Matjacic Z., Hunt K.J., "Paraplegic Standing Supported by 

FES-Controlled Ankle Stiffness", IEEE Transactions on Neural Systems and 

Rehabilitation Engineering, Vol. 10, December 2002, No. 4, pp. 239-248. 

[21] Holderbaum W., Hunt K.J., Gollee H., "Hcx) Robust Control Design for 

Unsupported Paraplegic Standing: Experimental Evaluation", Control 

Engineering Practice, Vol. 10, , 2002, pp. 1211-1222. 

[22] Hunt K.J., Gollee H., Jaime R.P., "Control of Paraplegic Ankle Joint 

Stiffness Using FES while Standing", Medical Engineering and Physics, Vol. 

23, 2001, pp. 541-555. 

[23] Van der Spek J.H., Veltink P.H., Hermens H.J., Koopman B.F.J., Boom 

H.B.K., "Static and Dynamic Evaluation of the Influence of Supplementary 

Hip-Joint Stiffness on Crutch-Supported Paraplegic Stance", IEEE 

Transactions on Neural Systems and Rehabilitation Engineering, Vol. 11, 

No. 4, December 2003, pp. 452-462. 

[24] Van der Spek J.H., Veltink P.H., Hermens H.J., Koopman B.F.J., Boom 

H.B.K., "A Model Based Approach to Stabilizing Crutch Supported 

Paraplegic Standing by Artificial Hip Joint Stiffness", IEEE Transactions on 

Neural Systems and Rehabilitation Engineering, Vol. 11, No. 4, December 

2003, pp. 443-451. 

[25] Kuo A.D., "An Optimal Control Model for Analyzing Human Postural 

Balance", IEEE Transaction of Biomedical Engineering, Vol. 42, Jan. 1995, 

pp. 87-101. 

[26] Don B. Chaffin, Gunnar B.J. Andersson, Bernard J. Martin, 

"Occupational Biomechanics", Third edition, John Wiley & Sons Inc., 1999. 

[27] Onyshko S., Winter D.A., "A Mathematical Model for the Dynamics of 

Human Locomotion", Journal of Biomechanics, Vol. 13, 1980, pp. 361-368. 

[28] Rapoport S., Mizrahi J., Kimmel E., Verbitsky O., Isakov E., "Constant 

and Variable Stiffness and Damping of the Leg Joints in Human Hopping", 

Journal of Biomechanical Engineering, Vol. 125, 2003, pp. 507-514. 

[29] Ferris D.P., Bohra Z.A., Lukos J.R., Kinnaird C.R., "Neuromechanical 

Adaptation to Hopping with an Elastic Ankle-Foot Orthosis", Journal of 

Applied Physiology, Vol. 100, 2006, pp. 163-170. 

[30] Farley C.T., Morgenroth D.C., "Mechanisms for Leg Stiffness 

Adjustment During Locomotion", North American Congress of 

Biomechanics, NACOB 1998. 

[31] A. Winter, International Society of Biomechanics, Biomedical Data 

Resources, Gait Data. Available: http : //www . isbweb . org/data/ 

[32] A.L. Hof, CGA Normative Gait Database, Limb Fitting Center, 

Laboratory of Human Movement Analysis, Department of Rehabilitation, 

University Hospital Groningen, Available: 

http://guardian.curtin.edu.au/cga/data . 

[33] "Optimization Toolbox for use with Matlab", the Math Works, Inc., 

User's Guide, Version 3, Available: www.mathworks.com . 



PROC. CAIRO INTERNATIONAL BIOMEDICAL ENGINEERING CONFERENCE 2006© 



